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SUMMARY 


Stokes’ stream function is studied for subsonic or supersonic flow 
past axisymmetric bodies of sma.ll slope . The first-order equation is 
found to be nonlinear. It can be linearized if one seeks only the formal 
order of accuracy of the slender-body approximation. In any case, serious 
loss of accuracy results from imposing the condition of tangent flow at 
the body on the mass flux rather than on the velocity. In a second 
approximation, neglect of the nonlinearity leads to a false solution even 
in the slender-body expansion. 


INTRODUCTION 


Stokes ’ stream function would appear to offer attractive possibilities 
for treating subsonic or supersonic flow past bodies of revolution of 
small slope. The condition of tangent flow at the body surface Is simpler 
for the stream function than for the velocity potential, because a condi- 
tion is imposed on the dependent variable itself rather than on Its 
derivative. Use of the stream function also permits streamlines away 
from the body to be determined easily. For these reasons, various investi- 
gators have adopted Stokes ’ stream function in linearized theory (refs . 1 
to 3 )- 

Another advantage might be anticipated if one proceeds to calculate _ 
higher approximations by iteration on the first, and considers practical 
shapes having stagnation points (rather than, say, cusped noses ) . Tien 
it has been found that the failure of the first approximation in the 
vicinity of the nose may In some cases be compounded so that the second 
approx 1 .ina.tion Is incorrect everywhere, and the danger of this is greater 
with the velocity potential 9 than with the stream function ijr . Thus 
in the analogous problem of plane flow, straightforward Iteration on the 
subsonic thin-airfoil solution for a sharp -edged profile (such as a 
biconvex section) yields a false second-order result even at mid-chord 
if one works with cp, but correct results (except very near the edges) 
if one uses \Jr (ref. 4 ). Both methods fail, however, for round edges. 

One can then abandon the thin-airfoil approximation, and by us ing exact 
conformal mappings find the proper second approximation for a round- 
edged airfoil (which must still, however, be corrected near the edge. 
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where it is not uniformly valid). Thus Eantzsche and Wendt treated the 
ellipse and Joukowski section by laborious computation, using the stream 
function (ref. 5)* Recently a second-order compressibility rule has been 
discovered that reduces any subsonic airfoil problem to an incompressible 
one; and since these difficulties disappear at zero Mach number, they are 
no longer an obstacle in plane flow (ref. 4). They would, however, 
reappear if one proceeded to a third approximation. 


For bodies of revolution, these matters are less thoroughly explored. 
The only relevant published solution is that of Schmieden and Kawalki who, 
working with Stokes * stream function, obtained a second approximation 
for subsonic flow past a slender unyawed .ellipsoid of revolution, and 
gave also the asymptotic expansion for Vanishing thickness that may be 
termed the second-order slender-body solution (ref. 6). They used con- 
formal mapping, so that analogy with Hantzsche and Wendt’s treatment of 
the ellipse would suggest that their treatment is valid. However, axisym- 
metric bodies induce weaker disturbances : than do their plane counterparts, 
so that stagnation points play a more dominant role, and questions of 
validity are accordingly more delicate. 


This is illustrated by the fact that, in contrast to the corresponding 
plane case, the second-order slender-body solution based upon cp fails 
at a round nose even in incompressible flow. Thus for the paraboloid 
whose radius is r = S\fx~, expansion of the exact incompressible solution 
gives to second order near the body 
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whereas straightforward iteration on the slender-body theory of Heaslet 
and Lomax (ref. 7) misses the last tern shown above. As in the plane 
problem for compressible flow, the missed term is an eigensoluticn 
associated with the blunt nose, being in fact the slender-body represen- 


tation of a point source at the origin, of potential 


i 5 4 (x 2 + r 2 )" 1 / 2 . 
16 


It induces no radial velocity on the axis, where the tangency condition 
is imposed, and so is overlooked if one works with 9. It does, however, 
Induce a mass flux, so that one might hope to catch it by working with 
t, and this expectation is in fact realized. One might therefore be 
encouraged to attack the corresponding cdmpressible problem using . 


For these reasons, the application of Stokes' stream function to 
compressible small -disturbance theory has been studied here in some 
detail. The results are surprising. Far; from offering the expected 
advantages , the stream function is found to “suffer from grave shortcomings . 
The proper first-order equation of motion is actually nonlinear. The 
nonlinearity can be disregarded in the first approximation if one seeks 
only the formal accuracy, as indicated by order estimates, of the 
asymptotic slender-body theory; but the numerical accuracy is in any case 
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far inferior to that obtained by using cp. In the second approximation, 
neglect of the nonlinearity is disastrous, so that Schmieden and Kawalki * s 
solution for the ellipsoid is in fact incorrect. 


SMALL-DISTUHBANCE EQUATIONS 


The continuity equation for axisymmetric compressible flow is satis- 
fied by introducing Stokes * stream function ¥, according to which the 
axial and radial velocity components in cylindrical coordinates ar e given 
by 


= fe!r ' 
P r 

= _B» 

P r 


CD 


where p is the density and its value in the free stream. Uniform 

inviscid flow past a body of small slope is irrotational to at least 
third order, and to that approximation the equation of motion is 


(c 2 - tF^xx + (c 2 - V 2 )frr - 2UV¥xr - c 2 = 0 (2a) 

where the speed of sound c is related to its free-stream value and 

the flow speed by 


c 2 + Z|i (ye + V 2 } = ^2 + Zfl ^2 


(2b) 


Consider now perturbations about a uniform flow, introducing pertur- 
bation quantities (denoted by lower-case symbols) by setting 


¥ = vj^ r 2 + *) 
U = Uco(l + u) 

V = U OT v 


(3) 


In terms of these perturbation quantities the equation of motion (2) 
becomes 

1 - M 2 (2u + u 2 + v 2 ) 


*xx + Tr- 


ier 


-jfc’ 

r 


) 


= f^[(l + 2u + u 2 )^ + V s (1 + T|r rr ) + 2(1 + u)vtxr] 
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The flow is isentropic to at least second order, and to that approxi- 
mation the density ratio appearing in equations (l) is given by 


ftp 

P 



1 _ Z_i + u 2 + v 2 ) 


1 

7-1 


( 5 ) 


Expanding the right-hand side for small u and v, and combining with 
equations (l) and (3) leads to the following expressions for the pertur- 
bation velocity components (referred to free-stream speed) in terms of 
first derivatives of the perturbation stream function: 
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(6b) 


( 7 ) 


Order of Perturbation Quantities 


It is well known that the streamwise velocity perturbation induced 
by a body of revolution is, except near porners , of smaller order than 
the radial velocity. For a body of thickness ratio r the radial com- 
ponent v is of order t as indicated by the tangency condition, but 
u is only of order r 2 Zn t. It follows: that the quantities (i|r r /r ) 
and (tx/ r ) 2 bracketed in the first term bf equation (6a) are o(r 2 Zn t) 
and 0 (t 2 ). Both are therefore of first order (even though one is linear 
and the other quadratic in perturbation quantities) because for slender 
bodies, terms of different order differ by factors of t 2 . Likewise, the 
terms in the second bracket are 0(t 4 , r 4 Zn r, T 4 Zn 2 T) and hence all of 
second order. Near corners u is 0 (t) because the flow is locally plane, 
and then the linear terms in equations (6) are of first order, the quad- 
ratic terms of second order, and the cubic and quartic terms are simply 
negligible . 

The orders of \|r and its derivatives near the body can be determined 
from equations (6) together with the facts that r is 0(r) and that 
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integration or differentiation with respect to x does not affect the 
order of a term. Aside from powers of In t, which are 0(l) for all 
practical purposes, the result is 


u = 0 (t 2 ) 


v = 0(t) 



Ofr 2 ) 

fxx 

= 0(t 2 ) * 

0(t) 

txr 

= 0(t 3 ) 

0(t 2 ) 

trr 

= 0(t 2 ) . 


( 8 ) 


First-order and Linearized Equations 


Substitut ing the expressions (6) for u and v into equation (4) and 
retaining only leading terms (of order t 2 ) yields the first-order equation 
of motion: 

(1 - + irr - ^jr = M 2 (^J (9) 

An unexpected result is that the first-order equation is nonlinear. 

All previous investigators (refs. 1 to 3, and 6) have based their 
first approximation upon the linearized equation obtained by dropping the 
right-hand side of equation (9): 

(1 - M 2 )*xx + t rr " T 1 = 0 (10) 

and have also linearized the connection between and u by omitting the 
quadratic first-order term in equation (6a). This procedure actually 
yields results correct to first order, as can be verified by introducing 
a modified perturbation stream function ijr according to 

t = ? - \ M 2 ? x 2 Zn ** (11) 

The order est ima tes of equations (8) apply to ijr as well as i|f. Hence 
substitution shows that the first-order problem for f implies, to first 
order, the linearized problem for \jr. However, the relation (6a) con- 
necting u and 1 Jr must be linearized as well as the equation of motion 
(10); to linearize either alone would lead to incorrect results. This 
suggests that the linearization may be somewhat of a coincidence, which 
cannot be relied upon to have its counterpart in higher approximations . 

It should be noted that, aside from enormous mathematical simplifi- 
cation, linearization leads to a problem satisfying the similitude of 
Gothert * s rule (ref. 8), which is not true of the first-order problem 
based upon equation (9). 
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Two theories that are asymptotically of equal order may differ con- 
siderably in numerical accuracy. Whether the nonlinear first-order 
equation or the linearized equation yields greater accuracy, and how they 
compare with the more conventional treatment using cp are questions to 
be considered in a subsequent example . — 


Second-order Equation 


Retaining terms of Ofr 4 ) in equation (4) yields the second-order 
equation of motion which, written in a form suited to iteration on the 
first-order solution, is i 


(1 - M 2 )txx + trr - --r - = M 2 


( 12 ) 

Here the right-hand side has been simplified by using the fact that the 
left-hand side is to be equated to zero in the first approximation. 

The counterpart of linearizing the first-order equation is the 
retention here only of quadratic terms, which yields, in a form suited 
to iteration on the linearized solution — 



(1 


M 2 )fex + 7rr 


r 


2(n + 



5 



fcr + 



(13) 

and this Is the basis of Schmieden and Kawalki’s second-order solution 
for subsonic flow past an ellipsoid of revolution (ref. 6). Corre- 
spondingly, they retain only quadratic terms in the expressions (6) for 
velocity components . 


FIRST APPROXIMATIONS FOR CONE IN SUPERSONIC FLOW 


It happens that the first-order equation (9), though nonlinear, can 
be integrated in closed form for the standard example of supersonic flow 
past an unyawed circular cone. This example can, therefore, serve as a 
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test case for assessing the relative accuracies of the varicrus possible 
first approximations . The proper choice of the relation between pertur- 
bation velocities and pressure coefficient is a possibly controversial 
question that will be avoided by restricting attention to the streamwise 
velocity increment on the cone. 


First-order Solution 


For flow past a cone, the stream function is homogeneous of order 2 
in the space coordinates, so that 


t(x,r) = r2f(t) 


(l4a) 


where 


t 

B 2 


= Br 

x 

= M 2 - 1 


(14b) 


Hence the first-order equation of motion (9) becomes the nonlinear 
equation for f(t): 


(1 - t 2 )f" + f » = ^ t 2 ! * 2 

t B 


(15) 


This is Bernoulli's equation for f which is a linear equation for 
l/f*, with solution 


f*(t) 


B 2 Jl-t g 
vF t 3 (k +sech -1 t) 


(16) 


The integration can be carried out in terms of exponential integrals. 
Imposing the condition that the disturbances vanish at the Mach cone, 
f(l) = 0, leads to 


f(t) 


_ 1 Bf C 
~ 4 M 2 X 


22n(^l 


sech""^ 


- e 


- 2 k 


Ei(2k + 2 sech“ X t) - Ei(2k) 


2 k 

e 


Ei( 


-2k - 2 sech 1 t) 


Ei(-2k) 


( 17 a) 
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where 


00 x 

ei(-x) = - yv ds * e!(x). = y> 


ds , x > 0 


( 17 b) 


First-order Solution With Tangency on Mass Flux 


The condition of tangent flow at the surface can now he imposed on 
either the mass flux or the velocity. In the full theory these are of 
course equivalent, hut in the approximate theory they yield results dif- 
fering in higher-order terms . Zero mass flux through the surface is 
assured hy requiring the stream function to vanish there. For a cone of 
semi vertex angle e this yields the transcendental equation 

f (Be) = - \ "V (17C) 

for determining k. Then the surface increment in streamwise velocity 
is given hy 


us 


^Js 1 _ 

Ux> B 2 


-1 + tf* 


ri§> <«■>* 


Jt=Be 


(17a) 


The transcendental equation (17c) has been solved numerically for a 
cone of 10° semi vertex angle. The resulting variation of u with free- 
stream Mach number is compared in sketch (a) with the exact results from 
Kopal’s table (ref . 9 ). 



A curious feature of the solution 
is that beyond a certain Mach number 
the constant k becomes negative. 

If ;k is slightly negative the factor 
(k + s^clf-H;) in equation (l6) will 
vanish on the cone t = sech(-k) which 
lies within the flow field, so that 
f 1 ,is singular there. This would 
mean that the radial velocity v was 
6 infinite in the flow field, which is 
intolerable . The solution must 
accordingly be regarded as having 
broken down when k becomes zero. 

There is no apparent physical interpretation of this limit . It arises at 
a lower Mach number than the more familiar cutoff of linearized theory at 
Be = 1, which corresponds to the free-stream Mach cone’s having been forced 
down onto the surface. Thus, for the 10° cone, k vanishes at Be = 0.699 
or M = 3-92, as indicated by the cutoff in sketch (a). 


I 


M 


Sketch (a) 
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First-Order Solution With Tangency on v 


Instead of requiring the stream function to vanish, one can impose 
the tangen cy condition on the velocity at the surface. In the approxi- 
mation us uall y adopted in linearized theory, the small increment in 
streamwise velocity is neglected, so that the condition is simply that 
v equal the slope of the hody. For the cone, this determines the 
constant k according to 


k _ Jl-B S e 2 

k ’ M 2 ^ 

Then on the surface 

1 + 2f (Be) + | M 2 ^ 

(18b) 

This is compared in sketch (b) with 
the exact solution for a 10° cone. 
The approximate solution is termi- 
nated at Be = 1, which occurs at 
M = 5-76. 


- sech“ 1 Be (l8a) 



M 

Sketch (b) 



Comparison With Other Approximations 


The linearized stream function for a cone has been found by Moore 
(ref. 3) by solving equation (10) and requiring ¥ to vanish on the 
surface. The streamwise velocity increment on the surface is 


us = -e 2 


sech" :L T 

*7 l-T 2 -T 2 sech”' L T 


( 19 ) 


where T = Be. Sketch (c) compares 
this result with the exact solution. 
Alternatively, one can impose the 
tangency condition on the velocity, 
which gives (with the approximate 
condition v = e discussed 
previously) 


Ug = -e 2 


sech" 1 ! 

n/i-T 2 


( 20 ) 



M 

Sketch (c) 
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The comparison with exact theory is 
given in sketch (d). 

Instead of working with the stream 
function, one can use the perturbation 
velocity potential 9, for which the 
linearized or first-order equation is 


(l - M 2 )cf>} QC + % r + - 0 


M 

Sketch (d) 


Then the approximate tangency condition v = e leads again to equa- 
tion (20). Alternatively, one can require zero mass flux through the 
surface. This condition has been adopted for mathematical simplicity by 
Parker in his analysis of minimum drag bodies (ref. 10), using a cylindri- 
cal control surface . The result is then equation (19) • Thus it appears 
immaterial whether one uses the linearized equation of motion for cp 
or \[rj the result depends only upon the farm of tangency condition adopted. 


It can be verified that each of the approximations considered so far 
has the same asymptotic form for small T = Be, namely the slender-body 
solution, according to which i 


ug = -e 2 ln T + 0(e 4 2n 2 T) (21) 

This approximation is compared with the exact solution for a 10° cone in 
figure 1, together, with all the other first approximations discussed 
previously. The various approximations differ only in the subsequent 
terms of the asymptotic series, which are not correct in any case because 
they are affected by the neglected nonlinear terms . Despite occasional 
statements to the contrary, however, this does not rule out the possibility 
tha t one approximation is consistently better than another. 

The most striking result is the inaccuracy associated with imposing 
the tangency condition on mass flux, whether in the linearized equation 
for 9 or f (eq, (19) and sketch (c)) or the nonlinear first-order 
equation for i|r (eq. (I7d) and sketch (a)). The error is never less 
than twice that in simple slender-body theory, and rises rapidly with 
Mach number. 

On the other hand, imposing the tangency condition on velocity leads 
in every case to results more accurate than those of slender-body theory. 
Although the nonlinear solution is the most accurate at high Mach numbers 
(eq. (l8b) and ske'tch (b)), it falls behind_at lower speeds; and in any 
case it could not reasonably be extended ito other bodies of revolution. 
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There remains the linearized solution (using either cp or t) with 
tangency imposed on v (eq. (20) and sketch (d)) as the most accurate 
approx ima tion applicable to other shapes . Experience has shown that it 
consistently yields accuracy of the sort encountered in this example. 

Hence it seems unlikely that Stokes’ stream function offers any advantages 
in first-order compressible flow theory, except possibly for determining 
streamlines in subsonic flow, where the inaccuracy would probably be no 
greater than it is at a Mach number of n^ 2 in supersonic flow. 


HIGHER APPROXIMATIONS 


It was suggested in the introduction that the stream function might 
have the advantage in higher approximations of applying to blunt bodies 
for which the velocity potential fails . However, a practical iteration 
scheme could probably only be based upon the linearized approximation. 
The discovery that the true first-order equation of motion is act ual ly 
nonlinear casts some doubt upon the validity of such a scheme. This 
matter will be studied by considering again the special case of a cone, 
for which the correct second-order solution is known. 

Iteration on Linearized for Cane 


The second approximation for a cone in supersonic flow will be sought 
by proceeding in strict analogy with Schmieden and Kawalkihs treatment of 
an ellipsoid in subsonic flow (ref. 6). Thus the slender-body expansion 
is introduced only at the last stage, which eliminates the well-known 
difficulty with the distant boundary condition in slender-body theory. 

With the previous substitution (l4) for the perturbation stream 
function, the linearized equation (10) becomes 

(1 - t 2 )f 1 " + 3 _ f 1 » = 0 ( 22 ) 

t 

which is just equation ( 15 ) without its nonlinear term. The solution 
that vanishes at the Mach cone t = 1 is 


(23a) 

The condition that the full stream function vanish on the cone of semi- 
vertex angle e deter min es the constant A as 

A = — (23b) 

\l 1-T 2 - T^sech” 1 ! 


f, = “ =■ 


</l-t £ 


- sech 


S ) 


These are just Moore’s results (ref; 3)- 
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Substituting this solution into the right-hand side of equation (13) 
gives as the iteration equation for the second -order increment to f 


(1 - t 2 )f 2 " + 3 - ~ 2t - g f 2 * 


M 2 A 2 

B 2 


2(E - _ 1 ) 


sech -1 t 


l+t £ 


t^s/l-t 2 t4 J 


(24a) 


where 


N 


_ 7+1 M 2 _ y-j-1 "M 2 

2 B 2 2 M 2 -l 


(24b) 


The general solution vanishing at the Mach cone is 

M 2 A 2 




sech _1 t ) + 


4b! 2 


si 1-t 2 sech -1 t 


(2W + 73) 


1-t 2 

t 2 


+ (l - 2N)(sech“ 1 t ) 2 


(25a) 


At this point it is appropriate to introducejthe slender-body approxi- 
mation by expanding asymptotically for small t. Then requiring the 
stream function to vanish on the surface to second order detenu ines the 
constant C as ' : 


C = M 2 e 2 T 2 (^ln | : W 



(25b) 


The velocity components are related to the stream function by 
equations (6). Schmieden and Kawalki’s analysis implies (falsely) that 
linear and quadratic terms are respectively of first and second order, 
which means that the streamwise velocity increment is given by 


u 




M 2 

2B 2 


"2N-5 (**r\ Z . /V\ 2 

. B 2 ; ) \ r ) 


(26) 


On the surface of the cone, the linearized solution for f x gives, in 
the slender-body expansion ~ ~ 



T 2 in | + T^ln 2 | + | In T 



(27a) 



r 


-e 


1 + T^ln =■ + 


T 



(27b) 
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and the second-order increment f 2 gives 

+ . . . (27c) 

Note the apparent contradiction that this supposed second-order term 
actually starts off with a first-order term in e 2 . That contribution 
to u is, however, cancelled by the first-order nonlinear term in 
(i)fi x /r) 2 , which is just the coincidence discussed previously that permits 
the first-order problem to be linearized. Thus equation (26) gives as 
the supposed second approximation 


¥2. 


- = - i M 2 e 2 + M 2 e 2 T 2 j 
r 2 


| - NjZn 2 | - I l n| + N + l 
2 / T 2 T 


us = 1 - e^ln 

Be 


B 2 Zn 2 — - 
Be 


i±ae 2 - + + i±as 

2 Be k 


( 28 ) 


This result is unf ortunately not entirely correct. Broderick's 
solution of the same problem, using the velocity potential, gives 6 m 2 /t 
in place of 3M 2 /4 in the last term. His solution has been confirmed 
by independent analyses, and must be regarded as reliable. That the 
discrepancy is in a term proportional to M 2 is plausible since the 
linearized and nonlinear first-order problems differ by terms in it 2 . 

If one imposes the tangency condition on velocity rather than mass 
flux, the result is different but still incorrect . A discrepancy then 
appears also in the coefficient of e^ln 2/Be. 

Discussion .- The preceding example shows that although the first- 
order problem for i|r can coincidental .1 y be linearized, the second-order 
problem cannot be limited to quadratic terms and solved by iteration on 
the linearized solution. Three modifications, of successively greater 
complexity, suggest themselves as the correct procedure. 

The simplest possibility is that the true second approximation can 
be found merely by retaining some cubic or quartic terms in the right-hand 
side of the iteration equation ( 13 ) or in the expressions (6) for velocity 
components . This is known to be the case if one works instead with cp, 
where the iteration equation must include a cubic term, and quadratic 
terms appear in the expression for pressure. However, attempts to isolate 
the proper combination (if it exists) from the many possibilities in the 
preceding example have not succeeded. Secondly, it may be that the 
remainder of the second-order terms would emerge from the next step of 
the procedure, involving cubic terms. Thirdly, it might be necessary to 
work with the nonlinear iteration equation (12), which would introduce 
enormous mathematical difficulties . 

The problem mentioned in the introduction of second-order subsonic 
flow past a paraboloid has recently been solved by the author in an indi- 
rectway, and the solution obtained by using ijr in the preceding fashion 
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is f oimd again to be in error by a term in M 2 . Since a paraboloid is 
a limiting case of an ellipsoid, this confirms the suspicion that 
Schmieden and Kawalki 1 s solution (ref. 6) is incorrect in the second-order 
terms in M 2 • ■ 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 10, 1956 
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Figure 1 ,- Various first approximations for surface speed on cone in supersonic flow. 
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